Adjustment of the electric current in pulsar 
magnetospheres and origin of subpulse modulation 
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Abstract The subpulse modulation of pulsar radio 
emission goes to prove that the plasma flow in the open 
field line tube breaks into isolated narrow streams. I 
propose a model which attributes formation of streams 
to the process of the electric current adjustment in the 
magnctosphcre. A mismatch between the magneto- 
spheric current distribution and the current injected by 
the polar cap accelerator gives rise to reverse plasma 
flows in the magnetosphere. The reverse flow shields 
the electric field in the polar gap and thus shuts up 
the plasma production process. I assume that a circu- 
lating system of streams is formed such that the up- 
ward streams are produced in narrow gaps separated 
by downward streams. The electric drift is small in 
this model because the potential drop in narrow gaps 
is small. The gaps have to drift because by the time a 
downward stream reaches the star surface and shields 
the electric field, the corresponding gap has to shift. 
The transverse size of the streams is determined by the 
condition that the potential drop in the gaps is suffi- 
cient for the pair production. This yields the radius of 
the stream roughly 10% of the polar cap radius, which 
makes it possible to fit in the observed morphological 
features such as the "carousel" with 10-20 subbeams 
and the system of the core - two nested cone beams. 
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1 Introduction 

The nature of pulsar emission is still a subject of 
debate. The conventional wisdom is that the pul- 
sar activity is associated with relativistic plasma flows 
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along the open magnetic field lines. According to 
the commonly accepted picture, the polar cap cascade 
produces electron-po si tron p airs (se e , e.g.. reviews by 
Grenier and Hardin j |2006l) : lAronsl (|2007l) and refer- 



ences therein), which stream along the open magnetic 
field lines with relativistic velocities. The pulsar radio 
emission is believed to be generated in this outflow. 

Because of strong relativistic beaming, the basic 
morphological features of the pulsar radiation are dic- 
tated by the geometry of flow. It is generally accepted 
that the width of the mean pulse is determined by the 
width of the open field line tube in the emission region 
so that narrow pulses imply low emission altitude. The 
subpulse structure is commonly attributed to a systems 
of narrow streams of plasma, each stream radiating a 
subbeam tangentially to the magnetic field lines. These 
subbeams exhibit pretty regular behavior. It was found 
that p ulsars have attracto rs in the form of two nested 
cones ( RankinI 1983 . 1993 ). Within the cone, the sub- 
pulse behavior also ex hibits substant i al regularity such 



Weltevrede et al 



as st a tiona ry drifting ()Rankinlll98 
20061 120071 ). In a couple of C8LS6S, cL careful analysis 
of the data reveals a remarkable circular "carousel" 



of emitting subbeams (jDeshpande and RankinI 12001 



Mitra and RankinI 120081) . We are still lacking of a firm 



model unequivocally explaining why the flow in the 
open field line t ube breaks into isolated streams. 

According to Ruderman and Sutherland! ( 1975 ). the 
electric field is unable to extract ions from the surface 
of the neutron star therefore electron-positron pairs are 
produced in local gap discharges (sparks). Each spark 
gives rise to a narrow plasma stream responsible for 
a subpulse. The sparks are isolated one from another 
because the potential drop in the space between them 
remains below that necessary for the plasma produc- 
tion. This model remains the most popular ( s ee, e.g ., 



recent paper by Ivan Leeuwen and Timokhinl (|2012f ) ) . 
The basic problem with this model is that the sur- 



face binding energy was found to be not high enough 
to prevent a free supply of charges from the surface. 
Then the primary beam is slowly accelerating so that 
the gap height exceed s the r a dius of the polar cap 



the charge density ( Goldreich and Julian 19691) 



(Muslimov and Tsvganl Il992t iHibschman and Arons 



2001al lb); in this case the electric se p aratio n of sparks 



does not work. iGil and Melikidzd (|2002r ); iGil et al 



( 20061 ) propose to resolve this difficulty assuming that 
the surface magnetic field is significantly larger than 
th e dipole . 



Wrightl (|2003l) proposed an empirical model for sub- 



pulse modulation based on the assumption that in the 
open field line tube, particles stream in both directions 
between the polar cap and the outer magnetosphere 
and these flows are separated one from another forming 
a slowly circulating system of plasma columns. With 
a number of other assumptions, this model nicely de- 
scribes the observed subpulse behavior however, the 
model is still lacking of a solid physical background. 
Wright suggested that the reverse flows are produced 
by the pair cascades in outer gaps. It is not at all evi- 
dent that outer gaps exist and could produce pairs, es- 
pecially in slowly rotating pulsars, and in any case, it is 
not clear why the plasma flows from the outer and polar 
gaps should be separated. However, the principal idea 
that the formation of the subpulse system should be at- 
tributed to an intercommunication between the outer 
magnetosphere and the polar cap sounds reasonable. 
An important point is that one need not look for an ac- 
tive region in the outer magnetosphere in order to find 
a feedback mechanism. Generally the structure of the 
magnetosphere is not driven by plasma producing gaps, 
neither by the polar nor by the outer ones, but on the 
contrary, the global electrodynamics dictates the charge 
and current distributions in the magnetosphere so that 
a gap has to adjust itself to produce the necessary cur- 
rent. How could the polar cap maintain the global 
magnetospheric current is still not clear; the problem 
is widely discus s ed in the literature ( see recent papers 
bvlAronsl(l2007l):lBeloborodovi|2008l ); lLuo and Melrose! 
(|2008| ); lTimokhinl (|2010l )). In this paper, I suggest that 
the plasma flow breaks into separate narrow streams as 
a result of the electric current adjustment in the open 
field line tube. 



2 Electric current adjustment in pulsar 
magnetospheres 

Within the magnetosphere, the plasma shields the lon- 
gitudinal (along the magnetic field) component of the 
rotationally induced electric field providing at any point 
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This condition is a generalization of the standard con- 
dition of the quasi-ncutrality because in the corotat- 
ing frame, the electric field vanishes when the charge 
density is equal to that of the Goldreich- Julian. Any 
deviation of the local charge density from (TTJ) results 
in a longitudinal electric field, which redistributes the 
charges to establish p — pej ■ Therefore in the plasma 
filled magnetosphere, the deviation of the charge den- 
sity from that of Eq. (fTJ) remains negligibly small at any 
point; one can say that the plasma shields the electric 
field in the corotating frame. 

The electric current within the open field line tube 
is also not a free parameter. The current distribu- 
tion is dictated by the global structure of the magneto- 
sphere. In the magnetosphere, the magnetic field lines 
bend backwards with respect to the rotation direction 
in order to allow the plasma to stream along them at 
a velocity smaller than the speed of light even though 
the rotation velocity becomes supcrluminal beyond the 
light cylinder. The current should be distributed such 
that the necessary magnetic field is maintained. For- 
mally, the current along any magnetic field line is de- 
termined by the condition of the smooth transition of 
the flow along this field line through the lig ht cylin- 
der (|lngrahamlll973t IContopoulos et al.lll999l) . Within 
the plasma flow, the necessary current density is main- 
tained owing to a small difference in the average elec- 
tron and positron velocities. The current flows along 
the magnetic field lines and the same current must 
flow across the polar gap at the base of the open field 
line tube. However, the flow in the gap is charge sep- 
arated so that there is no additional degrees of free- 
dom to adjust the current. The gap is filled only with 
the relativistically moving primary particles with the 
density close to pgj/s. Therefore the polar gap mod- 
els predict that the current density across the gap is 
close to PgjC- Even though the current density in the 
self-consistent magnetosphere is of the order of Pgjc 
, it is not equal to this value; moreover the current 
density is not the same over the polar cap. For ex- 
ample in the axisymmetric magnetosphere, the cur- 
rent density is equal to pgjc at the axis of the pul- 
sar, decreases towards the periphery of the polar cap, 
changes sign near t he boundary and again becomes zero 
at th e boundary ( Contopoulos et al. 19991 Timokhin 
20061) . In the oblique case, the current has opposite 
signs in the upper and low er part of the polar cap (e.g., 
Bai and Spitkovskvl (|2010h ). 

The mismatch between the current in the gap and 
the current flowing along the same magnetic field line 
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in the magneto-sphere means that "extra" charges are 
accumulated somewhere above the polar gap until the 
longitudinal electric field arises sufficient to drive back 
extra charges. Since the particle energy in the flow is 
small as compared with the total rotationally induced 
potential, the necessary excess of the charge density 
over the Goldreich- Julian density is very small. Thus, 
one sho uld expect that a reverse particle flow is easily 
formed (jLvubarskiil 119921 ); simulations of the plasma 
flow wit hin the open field line tube confirm this con- 
jecture ( Lvubarskv 2009h . The reverse particle flow 
fills the gap from above and shorts out the longitudi- 
nal component of the electric field provided the sign of 
the particle charges is the same as the Goldreich- Julian 
charge. In nearly aligned magnetospheres, this occurs 
at the periphery of the polar cap where the magneto- 
spheric current den sity is less than Pg.jc or even has the 
oppo site direction (|Contopoulos et al.lll999t iTimokhin 
20061) . Therefore when the reverse flow reaches the sur- 
face of the star, the particle acceleration is terminated 
and the pair production ceases. 

This could result in a non-steady regime such 
that the required current is maintained only on aver- 
age. Recently one -dimensional models w e re developed, 



both analytical dLevinson et al.l 120051 : iBeloborodov 
2008 ; ILuo and Melrosell2008j) and numerical (jTimokhin 
20101 ). for such a regime. These models assume that the 
gap height, h. is less than the polar cap radius and that 
the magnetospheric current is already fixed just above 
the gap. In this case one can imagine that different 
parts of the gap oscillate independently thus providing 
the necessary current distribution over the polar cap. If 
the gap height is larger than the polar cap radius, as it 
should be if the surface freely s upplies the charges (e.g 



Hibschman and Aronsl (l200rbl ): lHarding and Muslimov 

a5). 



|2002l ): lHarding et al.l l|2002l )), the situation should be 



much more complicated. Moreover, the formation time 
of the reverse flows, ~ O -1 , is much larger than the 
characteristic gap time and is generally different at dif- 
ferent magnetic field lines. In this case, one can hardly 
imagine that the necessary current distribution over the 
polar cap is achieved via the gap oscillations. 

Another possibility arises if one abandons the one- 
dimensional picture at all. Namely, the cascade could 
become highly inhomogeneous across the polar cap such 
that the required current distribution is achieved only 
by averaging over the small scale structure. One can 
speculate that the primary beam breaks into many sub- 
beams separated by streams of the reverse plasma flow. 
The primary subbeams are accelerated in the electric 
gaps and produce electron-positron plasma streams; 
the reverse flows are formed in the outer magneto- 
sphere via the charge and current adjustment. When 



a narrow reverse flow reaches the surface, the electric 
field is shielded at these field lines therefore the gap 
shifts so that the plasma is now produced at neighbor- 
ing field lines. This implies a slow drift of the entire 
system of the upward and downward moving plasma 
streams (Fig.l). Recall that in the axisymmetric mag- 
netos phere, the current density is equal to pg j c at the 
axis (jContopoulos et alJll999t ITimokhin! 120061 ) so that 
the reverse plasma flows arise only at the periphery 
of the polar tube. Therefore in pulsars with small in- 
clination angles, the flow breaks into a steady central 
stream surrounded by upward and downward streams 
slowly circulating around the magnetic axis. The pri- 
mary subbeams carry the current density pgjc and so 
do the plasma streams produced above the gaps. The 
reverse plasma flows could carry any current of any sign 
provided the pair density within them is at least a few 
times larger than the Goldreich- Julian density. If these 
streams are narrow enough, such a structure could pro- 
vide any current distribution averaged over the scale 
larger than the stream width but smaller than the width 
of the open field line tube. 

Even th o ugh t his structure resembles that proposed 
by Wright! ( 2003 ). the physics is different. Wright as- 
sumed that the potential drop in the gaps remains the 
same as in the Rudcrman-Sutherland model. However, 
the reverse plasma flows shield the electric field there- 
fore the potential drop in the gaps should be much 
smaller. In this case, the E x B drift rate is very low 
and could be neglected; the pattern circulates because 
the gaps could not remain at rest: they have to shift 
before the corresponding reverse flows come to the sur- 
face. Therefore the subpulse drift is determined not 
by the E x B drift rate but by the time necessary for 
the reverse flow to form and to reach the surface. The 
direction of the subpulse drift may be arbitrary. 

An important point is that a low potential drop in 
the gap does not generally preclu de the pair forma- 
tion. It was realized some time ago (IKundt and Schaai 



1993; iSturner et al.l Il995t iLuol 1996; Z hang and Qiao 



19961 ) that the inverse Compton scattering of the ther- 
mal radiation from the star surface provides pair- 
producing gamma-rays at the energy of the primary 
electrons much smaller than a few TeV required by 
the classical models. Depending on the structure 
of the magnetic field, the energy 1-10 GeV is suffi- 
cient in order to produce, in the Klein-Nishina regime, 
gamma photons, which are ea si ly converted into pairs 
( Hibschman and Aronsl2001allb : Harding and Muslimov^ 



2002[lHarding et al.ll2002l : lArendt and EilekbooilMedin and Lai| 



2010p. In this case, the pair yield is significantly smaller 
than in the classical models however, the plasma den- 
sity larger than the Goldreich- Julian density (multi- 
plicity larger than unity) could be achieved at typical 
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Fig. 1 A view of the plasma flows in the pulsar magnetosphere in the frame of the star. The flow breaks into a system 
of narrow streams. Each outward stream is produced in a corresponding potential gap where the primary particles are 
accelerated and give rise to the pair cascade. When the stream reaches the outer magnetosphere, a fraction of particles is 
redirected back because of the electric current mismatch between the polar gap and the magnetosphere. When the reverse 
plasma flow reaches the surface of the star, the electric field is shielded therefore the gap has to shift, which results in the 
drift of the whole pattern (to the left in the picture). All the particles move along the magnetic field lines (strictly vertical 
in the picture) however, the plasma streams are bent because the polar gaps (particle sources) drift so that the plasma 
produced at some bundle of the magnetic field lines reaches large altitudes when the plasma producing cascade is already 
shifted to the neighboring magnetic lines. 
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pulsar paramet ers if the magnetic field near th e surface 
is not a dipolc ( Hibschman and Aronsl 2001b ). 



3 Structure of the multi-beam polar cap 
accelerator 

In the picture outlined above, the polar cap region is or- 
ganized into an assembly of narrow gaps producing the 
outwards plasma streams. The gaps are separated by 
reverse plasma flows. Since the reverse flows shield the 
electric field, the potential drop in the gaps is small, the 
narrower the gap the smaller the potential drop. The 
condition that the gaps produce the electron-positron 
plasma places a lower limit on the gap size. In this sec- 
tion I demonstrate that the radius of the gaps could be 
as small as 10% of the polar cap radius; then the num- 
ber of steams is compatible with the observed picture 
of drifting subpulses. 

Even though a large potential drop is not neces- 
sary for pair production, the acceleration process could 



braking force (Kardashevet al. 1984; 


Xia et al.l 


Il985l 


Daughertv and Harding 1989: 


Chanel 199,4 


Sturnerl 


19951). The strongest is the force due to the resonance 



cyclotron scattering. Let the surface of the star radi- 
ate as a black body with the temperature T and the 
particle move perpendicularly to the surface with the 
Lorentz f actor 7. Then the resonance braking force is 
found as jDermerlll990T ) 



F 



c 2 7 



In 



(2) 



where is the cyclotron frequency, r e the classical 
electron radius. The maximal force is 



Q 4g r e cJBfc T = o 057Sl2T 2 meC 2 cm -i. (3-, 



where B = 10 12 £?i2 G is the surface magnetic field, 
Tq = T/10 6 K. It is achieved at 7 = 0.7haj B /kT = 
100i?i2Tg~ 1 when the cyclotron frequency falls into the 
Lorentz shifted thermal peak. The primary particles 
could acquire larger energies if the electric field in the 
gap exceeds F max /e. 

When the primary particles are extracted from the 
surface of the star, they are accelerated at first as in a 
plane-parallel slab. The electric field in this zone arises 
du e to the spac e charge limited effect; it is estimated 
as (|Miche]lll97i 



eE = ^/8TTep GJ m e c 2 = 0.3B{^ 'P~ 1/2 m e c 2 cm" 1 . 



(4) 



One sees that at typical parameters, the electric force 
exceeds the braking force ©. This field is extended 



up to the altitude comparable with the radius of the 
primary beam, a, beyond which it decreases exponen- 
tially. In classical models, the radius of the primary 
beam is equal to the polar cap radius. In our model, 
the primary beam is split into narrow subbeams sep- 
arated by reverse plasma flows therefore in our case, 
a is the radius of a subbeam. The field of the space 
charge accelerates particles up to the Lorentz factor 
7 - eEa/m E c 2 = 3005^ p-^^a/lO 3 cm). 

Above this zone, the particle is accelerated because, 
due to the general relativistic frame dragging, a mis- 
match arises between the charge density of the primary 
beam and the local Gqldreich- Julian charge density 



( Muslimov and Tsvgan 
1997|), 



1992t iMuslimov and Harding 



PGJ 



2GI 



O.5/45 



c 2 (i?» + zf 



(5) 



Here i?* = 10 6 i?*6 cm is the radius of the star, / = 
10 45 /45 g-cm 2 the moment of inertia of the star, z the 
altitude above the surface (I assume that z <C i?*). 
Within a narrow primary subbeam surrounded by con- 
ducting reverse plasma flows, the electric field is nearly 
transverse so that the electric potential is estimated as 
<f> = 2irApa 2 . Then the accelerating electric field is 



dz 



a Pgj- 



(6) 



Numerically the accelerating force could be presented 
as 

eE = Q^-(—\ m e c 2 cm- 1 . (7) 
P 2 \a J 



Here 



a 



(8) 



is the polar cap radius found for the vacuum dipolc 
field. Note that simulations of the pulsar magne 



19991 iTimokhinl 12006 



tospheres (iContopoulos et al . „ 

Bai and Spitkovskvl 2010 ) show that the real radii of 
the polar caps arc approximately 10% larger than ao- 
One sees that the accelerating force exceeds the brak- 
ing force ([3]) if a > 0.lPT 6 ao- Then the particles are 
accelerated up to large Lorentz factors sufficient for the 
pair production. 

Due to the frame dragging effect, the primary parti- 
cles could acquire Lorentz factors 
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eER* 



m e c 



= 6 • 10 4 



B 



12 



P 2 \0.la o 



(9) 



For small a, these electrons could not produce pairs via 
curvature radiation however, the inverse Compton scat- 
tering of the thermal radiation from the surface of the 
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star could initiate the pair cascade. The one-photon 
pair production occurs at e > 1000m e c 2 . The reso- 
nance scattering provides photons with the energy e = 
jhajs = 20074_Bi2?7i. e c 2 . At 7 > 10 4 , the non-resonant 
scattering in the Klein-Nishina regime produces more 
conv ertible photons than the r esonance scattering 
docs (|Hibschman and Aronsl 12001a ). Namely, the non 



reson ance braking force is (jBlumenthal and Gould 
1970h 



nNR 



71" 

y 



kT \ 2 a 2 2 / 4 7 fcT 

—m e c In - - 2 



m e c 



)2T|(ln^ 
1500 



— 2 I m P c z cm" 



(10) 



where a is the fine structure constant, Ac = h/m e c 
the Compton wavelength. In this regime, the en- 
ergy of the scattered photons is close to the electron 
energy therefore when the electron with the Lorentz 
factor 7 ~ fewxlO 4 passes the distance ~ it 
produces ~ F NB -R*/'ym e c 2 ~ 1 photon, which in 
turn gives rise to a few pairs with a smal l er ene rgy 
(for more deta i ls, see lllibschman and Aronsl ( 2001al lbl) ; 
Harding et al.l (|2002j)). One sees that even if the width 
of the primary subbeams is roughly ten times less than 
the width of the polar cap, the primary electrons pro- 
duce enough pairs to maintain the charge and current 
densities required by the pulsar magnctosphcrc. 

4 Radiation pattern and drift rate. 

Now let us come round to the overall picture. As a re- 
sult of the electric current mismatch between the polar 
cap accelerator and the magnetosphere, the flow in the 
open field line tube breaks into a system of the out- 
ward streams separated by the reverse plasma streams. 
The reverse streams reach the surface and short out 
the electric field forming plasma columns connecting 
the star with the outer magnetosphere. The electric 
gaps are formed between the reverse flows; the primary 
subbeams are accelerated in the gaps and produce the 
outward streams of the secondary plasma. The current 
density in the outward streams remains close to the 
Goldrcich- Julian value, as is dictated by the properties 
of the polar cap accelerator. The "extra" current flows 
in the downward streams. The narrower the streams, 
the closer the current distribution to that demanded by 
the magnetosphere therefore one can assume that the 
polar cap is populated as densely as possible with the 
streams. According to the above estimates, the radius 
of the subbeam could be as small as O.lao. The radius 
of the reverse flows should be approximately the same 



therefore the distance between the centers of adjacent 
subbeams could be taken as 0.4ao- 

In this model, the potential drop is small therefore 
the E x B drift could be neglected. The subbeams drift 
because they could not remain at rest as the electric 
field is shielded when the reverse plasma flow reaches 
the surface of the star. An important point is that in 
aligned magnetospheres . the current is equal to pgjc at 
the axis (|Timokhirj2006l ) therefore in pulsars with small 
inclination angles, the reverse flow is not formed near 
the axis of the polar cap so that the central subbeam 
could remain at rest. In this case other subbeams could 
circulate around the center of the polar cap. Taking 
into account the above estimate of the minimal distance 
between the subbeams, one can place the streams on the 
polar cap such that they form up to two concentric rings 
and a stream in the center of the cap. This is compat- 
ible with the basic observed radiation pattern: an ax- 
ial (core) beam surrounded by two nested cone beams. 
The number of streams in the outer ring could be esti- 
mated as 27rao/(0.4ao) = 16, which is compatible with 
the number of subbeams ( 20) fo u nd observationally 
by Deshpande and Rankinl ( 2001 '): Mitra and Rankinl 
(200S 



The circulation velocity is determined by the time 
necessary for the plasma stream reach the outer mag- 
netosphere, where the electric current mismatch pro- 
duces the reverse flow, and then by the time necessary 
for the reverse flow to reach the surface of the star. Let 
us estimate the time necessary to reach the outer mag- 
netosphere and come back. The particles move along 
the rotating magnetic field lines so that their velocity 
could be presented as 



v = fl x r + £B. 



(11) 



The constant £ is found from the condition that the 
velocity is close to the speed of light, which yields 



-B+x ± y'B 2 - B 2 x 2 
B 2 : 



(12) 



where x = Qr/c = 2itr/cP is the cylindrical radius in 
units of the light cylinder radius, P the pulsar period, 
B p and B$ the poloidal and azimuthal components of 
the magnetic field, respectively. The sign ± refers to 
outward and downward motion, respectively. Note that 
if fl and B p are positive, B^ is negative. The time 
necessary to reach a point at the cylindrical distance 
Xi from the axis or to come back from this point to the 
star is 



t±(a?i) = ± f 
Jo 



dx 
£±B X 



(13) 
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where the integral is calculated along the field line. The 
drift rate is determined by the condition that the gap 
shifts by its diameter, 2a, for the time t + + t—. It is 
customary to to characterize the drift rate by the time 
interval P 2 between two consecutive subpulses to pass 
the same longitude. In our model, P 2 = 2(i + +t-). 

The reverse flows are formed due to the electric cur- 
rent mismatch between the magnetosphere and the po- 
lar cap. Since the magnetospheric current distribution 
is determined by the conditions of s mooth transition 
of the flow through the l ight cylinder rtlngrahamlll973 ; 
Contopoulos et al. 19991) the reverse flows should be 
formed near the light cylinder. Note that £_ goes to 
zero at the light cylinder so that t^(xi) diverges as 
X\ — > 1 therefore the value of P 2 is determined by the 
details of the reverse flow formation. Just for the es- 
timate, I calculated t± for the flow along the bound- 
ary between the closed and open par ts of the magneto - 
sphcre in the axysimmetric model bv lTimokhinl ([20061 ) . 
The numerical data were kindly presented by Andrey 
Timokhin. I obtained P 2 = LIP for x x = 0.9, P 2 = 2P 
for x x = 0.95 and P 2 = 3.8P for x\ = 0.97. There- 
fore the typically observed P 2 ~ 2P could be easily 
obtained in the scope of this model. A significantly 
larger P 2 could be obtained only if the reverse flow is 
formed extremely close to the light cylinder. 



5 Conclusions 

In this paper, the subpulse modulation of the pul- 
sar radio emission is attributed to the current mis- 
match between the magnetosphere and the polar cap 
accelerator. The current distribution in the magne- 
tosphere is determined by the global electrodynamics. 
For example in the axisymmetric magnetosphere, the 
current density is equal to pgjc at the axis and de- 
creas es towards the peripher y of the open fi eld line 
tube |Contopoulos et all Il999t iTimokhinl 120061) . The 



problem is that the polar cap accelerator could inject 
only the current close to pajc- The "extra" charges 
should be expelled from the magnetosphere . which im 



plies the formation of reverse particle flows (jLvubarskv 
20091) . If the "demanded" current is less than the in- 



In this scenario, the plasma is generated in nar- 
row gaps scattered over the polar cap, which yields 
a system of narrow upward moving streams with the 
current density close to pgjc In the light cylinder 
zone, a fraction of this plasma is redirected back form- 
ing reverse plasma streams. When the reverse streams 
reach the surface of the star, the plasma production 
ceases at the corresponding bundles of magnetic field 
lines therefore the plasma producing gaps have to shift, 
which results in the slow circulation of the pattern. 
The current in the reverse streams is not restricted 
therefore the overall current distribution averaged over 
the narrow streams could easily match that demanded 
by the magnetosphere. In this model, the subpulse 
modulation arises naturally. The transverse size of 
the streams is limited from below by the condition 
that the plasma could be produced in the gaps. I 
demonstrated that the plasma production process is 
still possible if the radius of the streams is as small 
as 10% of the polar cap radius. In this case, 10 - 
20 subbeams could be placed along the circumference 
of the polar cap, which is co mpatible with the obser- 
vations of drifting subpulses ( Deshpande and Rankinl 
200lllMitra and Rankinll2008l ). Taking into account the 
estimated size of the streams, one sees that they could 
be arranged into a system of two concentric rings and 
a central stream, thus reproducing a generic radiation 
pattern: the core be am surrounded by two nested cones 
(|Rankin|[l983[|l993l) . 
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the results of his simulations of the pulsar magneto- 
sphere. This work was supported by the Israeli Science 
Foundation under the grant 737/07. 



jected one, the reverse flow brings back charges of the 
same sign as pgj- These charges fill the polar gap and 
shield the longitudinal electric field there thus arresting 



a spasmodic plasma production (Levinson et al. 


2005; 


Luo and Melrose!l2008; 


Timokhin 


I2OIOI). Here I suggest 



that the adjustment of the electric current between the 
magnetosphere and the polar cap accelerator occurs via 
the formation of narrow streams transferring different 
currents. 
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